Mon. Not. R. Astron. Soc. 000. ITU221 (2002) Printed 6 February 2008 (MN WT$t style file v2.2) 



PKS B1545-321: Bow shocks of a relativistic jet? 



V. Safouris 1 ' 2 *, R. Subrahmanyan 2 ' 3 , G. Bicknell 1 and L. Saripalli 2 ' 3 

1 Research School of Astronomy and Astrophysics, Mount Stromlo Observatory, Australian National University, 
Cotter Road, Weston ACT 2611, Australia 

2 Australia Telescope National Facility, CSIRO, Locked Bag 194, Narrabri, NSW 2390, Australia 
3 Raman Research Institute, C V Raman Avenue, Sadashivanagar, Bangalore 560080, India 



6 February 2008 



00 

o 
o 

(N 

in 



6 

CO 



> 

00 

m 
in 
o 

(N 
O 
oo 
O 



ABSTRACT 

Sensitive, high resolution images of the double-double radio galaxy 
PKS B1545— 321 reveal detailed structure, which we interpret in the light of pre- 
vious work on the interaction of restarted jets with pre-existing relict cocoons. 
We have also examined the spectral and polarization properties of the source, 
the color distribution in the optical host and the environment of this galaxy in 
order to understand its physical evolution. We propose that the restarted jets 
generate narrow bow shocks and that the inner lobes are a mixture of cocoon 
plasma reaccelerated at the bow shock and new jet material reaccelerated at 
the termination shock. The dynamics of the restarted jets implies that their hot 
spots advance at mildly relativistic speeds with external Mach numbers of at 
least 5. The existence of supersonic hot spot Mach numbers and bright inner 
lobes is the result of entrainment causing a reduction in the sound speed of the 
pre-existing cocoon. The interruption to jet activity in PKS B1545— 321 has been 
brief - lasting less than a few percent of the lifetime ~ (0.3 — 2) x 10 s yr of the 
giant radio source. The host galaxy is located at the boundary of a large scale 
filamentary structure, and shows blue patches in color distribution indicative of 
a recent merger, which may have triggered the Mpc- scale radio galaxy. 

Key words: galaxies: active - galaxies: individual (PKS B1545— 321) - galaxies: 
jets - radio continuum: galaxies 



1 INTRODUCTION 

PKS B1545— 321 is a remarkable example of a rare type 
of powerful radio galaxy in which two dou ble radio struc- 
tures share the same radio axis and core l|Saripalli et al.l 
120031 '): PKS B1545-321 is a double-double radio source. 
It was proposed that the inner 300 kpc double source, 
embedded within the Mpc-scale lobes of the giant radio 
galaxy, represents a new cycle of activity; that is, newly 
restarted jets have been 'caught in the act' of propagating 
through the remnant cocoon of a previous active phase. 
There are more than 10 such sources in the literature that 
have been identified as possible restarting radio sources 
on the basis of the detection of inner double structures 
recessed from the ends of diffuse outer lobes. In the case 
of PKS B1545— 321, the outer lobes extend almost all 
the way to the centre without any significant emission 
gap and the inner double is traced, embedded within the 
outer lobes, all the way from the central core to its bright 
ends. Therefore, PKS B1545— 321 is a good candidate for 
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a detailed study of the interaction of restarted jets with 
relict cocoons. Understanding the physics in such interac- 
tions is a step towards using the phenomenon of double- 
double radio sources as a probe of the physical nature of 
extended radio lobes and the causes and timescales asso- 
ciated with recurrence in jets in active galactic nuclei. 

The existence of double-double radio sources is evi- 
dence that the jet activity from the central engine of an 
active galactic nucleus (AGN) is not continuous over the 
life-time of a source: There appear to be interruptions to 
the flow that may be related to either a lack of fuel at 
the central engine, or instabilities in the accretion disk, 
or instabilities in the jet production mechanism. It has 
been recognized that unless the axis of the central en- 
gine changes dramatically, the newly restarted jets will 
initially propagate through the relict cocoons of past ac- 
tivity, rather than thermal ambient medium, as in the 
usual case of a single jet outburst. We expect the evolu- 
tion of the jets to be unusual in the cocoon environment, 
since the energy density in the lobes of a radio galaxy 
is dominated by relativistic magnetized plasma and the 
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cocoon medium is much lighter compared to the thermal 
interstellar and intergalactic media. 

Models for the development of jets in syn- 
chrotron cocoon s have been present e d in t he literature. 
IClarke fc Burnsl (|l99lh and Iciarkd (ll997T l carried out 
two- and three-dimensional numerical simulations of an 
under-dense, supersonic (Mach 6) restarting jet, and 
found that the new jet is heavy with respect to the co- 
coon of the first jet. As a consequence, the restarted jet 
propagates almost ballistically through the pre-existing 
cocoon, terminating in only a weak (and hence not very 
radio bright) shock. The supersonic advance of the new 
jet in the old cocoon also excites a weak bow shoc k imme- 
diately ahead of the new jet, which IClarkd (1 19971 ) suggest 
should be a visible feature in the radio l obe of a restart- 
ing ra dio so urce. Another feature of the lClarke fc Burnsl 
(| 19911 ) and IClarkd <|l997l) simulations is that there are 
no emission features that could be ascribed to emission 
from the cocoon of the restarted jet. This is presumably 
the result of the weak terminal shock combined with the 
rapid adiabatic expansion of the post hot spot flow. 

Whilst these simulations were informative and high- 
lighted aspects of restarting jets that were not appreci- 
ated at the time (in particular the rapid advance of the 
head of the new jet), they are not consistent with the 
observations of the Mpc-scale double-double radio galax- 
ies of which PKS B1545— 321 is an example. Bow shock 
features have not been realized in observations of restart- 
ing sources to date: The observed inner double sources 
simp ly appear as edge-bri g htened lobes in radio image s 
(e.g. ISaripalli et al.l (|2002h . ISchoenmakers et all i|2000l '), 
ISaikia et al.l (|2006T )). Moreover, the outer ends of the in- 
ner doubles have a high contrast with respect to the back- 
ground radio lobes indicating that the terminal shocks 
are at least of moderate strength. (We note however, an- 
ticipating observational results discussed below, that the 
inner hot spots are not as bright as one observes in nor- 
mal FR2 radio galaxies.) 

A different approach to the observations is to inter- 
pret the inner doubles as edge-brightened lobes, rather 
than near-ballistic jets and to invoke thermal densities 
inside the pre-existing cocoons at levels tha t are higher 
than predicted by simulations. iKaiser et al.l (2000) pro- 
posed that the entrainment of dense warm clouds of the 
intergalactic medium (IGM) could significantly contami- 
nate the synchrotron emitting cocoons, so that the devel- 
opment of the new jets thereafter would be more consis- 
tent with the observations: The interaction between the 
new jet and entrained gas would result in leading hot 
spots and new, filled inner lobes. Additionally, the long 
timescales for the entrainment and dispersion of the con- 
taminating material make this type of model particularly 
relevant to giant radio galaxies, in which inner doubles 
are frequently observed. However, leading bow shocks are 
expected in this model as well. 

Interpreting the inner structures as lobes may solve 
one problem, namely their brightness, but the problem 
of a missing bow shock remains. An additional possibil- 
ity is that the restarted jets may advance hypersonically 
with respect to the relict synchrotron lobes of the ra- 
dio galaxy. In this case, the leading bow shock may have 
a sufficiently small opening Mach angle that the emis- 



sion from both the bow shock and restarted jet would 
together form a narrow inner double structure, appear- 
ing as inner lobes without any distinct lead ing bow shock . 
This possibility was in fact considered by IClarkd (|1997h 
as a possible interpretation of the quasar 3C219. How- 
ever, the required Mach number was of order 1,000. In 
the case of PKS B1545— 321 the required Mach number 
is not nearly so extreme so that this idea is worth pursu- 
ing. Other variations on this theme include the possibility 
that changes in direction of the restarting jet would widen 
the working surface and lead to an apparent decrease in 
the bow shock angle with consequently less stringent re- 
quirements on the Mach number of the leading hot spot. 
Another possibility is that the bow shock is further in- 
fluenced by the cocoon matter distribution. We consider 
these possible theoretical interpretations in §[5] and espe- 
cially § 15.41 

In this paper we present new higher-resolution VLA 
radio continuum observations of the inner double source 
of PKS B1545— 321 at 22 and 6 cm wavelengths, and ex- 
amine the detailed structure of the inner lobes in the light 
of the models described above. Our modelling of the prop- 
agation and evolution of restarted jets in a relict cocoon 
results in estimates for the entrained thermal densities 
in the outer lobes and the external Mach number of the 
new jet. With about five synthesized beams across the in- 
ner double, our new VLA images of PKS B1545— 321 are 
the highest quality (sensitivity and resolution) images to 
date of any restarting source and provide an excellent op- 
portunity for such an investigation. Hydrodynamic simu- 
lations of restarting jets are being performed in order to 
provide an insight into the appearance of inner doubles 
in these scenarios; the results of these simulations will be 
presented elsewhere (Safouris et al., in preparation). 

We adopt a flat cosmology with parameters flo = 
0.3, Qa = 0.7 and a Hubble constant Ho = 71 km 
s _1 Mpc^ 1 . The host galaxy has an absolute magni- 
tude M B = -20-9 a nd is at a redshift of z = 0.1082 
l|Simpson et al.lll993h . At this redshift, l"= 1.95 kpc. 



2 RADIO CONTINUUM OBSERVATIONS 

Radio observations of the giant radio galaxy 
PKS B1545— 321 were previously carried out with 
the Australia Telescope Compact Array (ATCA) in 
2001 jSaripalli et al.ll2003l ): three separate 12 h Fourier- 
synthesis observations were done in three different array 
configurations at 22 and 12 cm wavelengths. 

We have made new observations of the source with 
the Very Large Array (VLA) in 2002, at 22 and 6 cm 
wavelengths, to improve the imaging sensitivity and res- 
olution for the inner source structure: the improved radio 
continuum images form the basis of this work. A journal 
of these observations is given in Table 1. The source, lo- 
cated at the southern declination of —32°, is visible to 
the VLA for a total of 6 hr in any single observation. 
We used the hybrid array configurations with the north- 
ern extended arm to give near circular beams for this 
southern source. To obtain the maximum resolution, we 
used the BnA array configurations at both observing fre- 
quencies. The source was also observed at 6 cm with the 
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scaled CnB configuration to obtain visibilities with spac- 
ings that match the 22 cm BnA data: multi-frequency 
observations using scaled arrays are important for stud- 
ies of spectral index and polarization as a function of 
wavelength. 

In the 22 and 6 cm bands, we used band widths of 
25 and 50 MHz centered at 1384 and 4910 MHz respec- 
tively. These bandwidths were divided into independent 
channels 3.125 and 12.5 MHz wide respectively, in order 
to avoid image degradation owing to bandwidth smear- 
ing. At 22 cm, we estimate that the reduction in peak 
response at the ends of the 9' outer lobes, as a result 
of bandwidth smearing, is less than 5%. At 6 cm wave- 
length, the reduction in peak response at the ends of the 
3'0 inner double is also less than 5%. 

The multi-channel continuum visibility data were 
calibrated using standard techniques in AIPS. The chan- 
nel data were averaged to make a pseudo-continuum 
'channel-0' data-set. Gains and phases in the 'channel-0' 
visibility data were calibrated using interspersed obser- 
vations of the nearby secondary calibrator J1522— 275. 
The flux scale was set to the Baars scale l)Baars et al.l 
using an observation of 3C286. The observation of 
3C286 was also used to calibrate the instrumental polar- 
ization. The 'channel-0' amplitude and phase calibration 
was copied to the multi-channel data and bandpass cal- 
ibration was determined using observations of the sec- 
ondary calibrator. 

The calibrated 6 cm BnA and CnB visibilities ob- 
tained using the VLA were concatenated. We constructed 
images from the 22 and 6 cm VLA visibilities sepa- 
rately using MIRIAD, adopting the Clarke algorithm 
(jClarkl Il980h to deconvolve the strong compact sources 
and the Steer-Dewdney-Ito algorithm (|Steer et al.lll984h 
for fainter extended source structure. Several rounds of 
phase self- calibration were iteratively performed. We 
combined the self-calibrated VLA 22 cm visibilities with 
self-calibrated archi val 22 cm ATCA visibilities (from 
ISaripalli et aT]|2003l ); to our knowledge this is the first 
time that VLA and ATCA datasets have been com- 
bined and imaged. We performed a linear mosaic of these 
22 cm VLA and ATCA visibilities using the MIRIAD rou- 
tine INVERT, which accounts for the differing primary 
beams of the two telescopes. After converting the ATCA 
linear and VLA circular polarization measurements to 
Stokes parameters, the routine images the two telescope 
data separately and then combines them in a linear mo- 
saic process, where pixels in the individual images are 
weighted to correct for the primary beam attenuation. 
Deconvolution of the mosaiced 22 cm image was done us- 
ing the maximum entropy algorithm implemented in the 
MIRIAD routine MOSMEM. Our technique for combin- 
ing and imaging the separate ATCA and VLA observa- 
tions is described in further detail in Appendix IA1 



2.1 22 cm radio continuum 

The combined VLA and ATCA 22 cm radio image is 
shown in Figure [1] The image was made with a beam of 
FWHM 3'.'6 x 2'.'9 and has a dynamic range of about 200. 
The rms noise in the image is a — 40 (iJy beam -1 and the 



lowest contour is drawn at 5a. The total angular extent of 
the radio source is 8.' 7 (corresponding to a projected lin- 
ear size of 1.02 Mpc). From the core, which is coincident 
with the nucleus of the host galaxy, the source is some- 
what more extended towards SE; the ratio of the extents 
towards SE and NW is 1.07. The low-surface-brightness 
lobes are edge brightened and the bridge between the 
ends is connected with no emission gap. There are no hot 
spots at the ends. The brightest feature in the image is 
the 2 .'9 (340 kpc) inner double structure, which appears 
to be fully embedded within the diffuse lobes of the giant 
radio galaxy. 

The lobes have a surface brightness that increases 
towards the ends: the center of the bridge has a sur- 
face brightness about 100 mjy arcmin -2 rising to about 
300 mjy arcmin -2 towards the two ends. Relatively 
bright rims and weak emission peaks are observed at 
the ends, which are just factors of two brighter than the 
surrounding lobes. These weak emission peaks, or warm 
spots, are well resolved by the 3'.'6 x 2 / .'9 beam and are 
likely the remnants of past hot spots: their presence sug- 
gests that the sound crossing times in the lobes are not 
high enough to erase the evidence of the powerful jets 
since they stopped feeding the ends of the outer lobes. 
The radio contours in Figure [T] which are spaced in loga- 
rithmic intervals, are tightly spaced at the ends indicating 
a sharp boundary; the relict lobes do not appear to have 
had the time to relax to equilibrium with the surrounding 
IGM. 

The total flux density of the outer double radio 
source is 1.7 Jy at 22 cm, which corresponds to a to- 
tal radio power of 5 x 10 25 W Hz -1 at this wavelength. 
Considering the absolute magnitude of the host galaxy, 
Mb = —20.9, this places the outer lobes of the Mpc-scale 
radio galaxy in the powerful FR-II radio galaxy regime 
(|Owenlll993l ). and almost an order of magnitude above 
the FR I/II break in radio power. This is consistent with 
the edge-brightened morphology of the giant radio lobes. 

The inner double structure is clearly seen in the 
22 cm radio image as a pair of relatively bright almost lin- 
ear features that are entirely contained within the outer 
lobes of the giant radio galaxy. In projection, the inner 
double and outer lobes have a common radio axis and 
core. The inner source has been traced in the 22 cm im- 
age all the way from the central core to its bright ends, 
which are well recessed from the outer lobe ends. A grad- 
ual curvature is observed over the entire length of the in- 
ner source, which is reflection symmetric about the core. 
An extrapolation of the inner double structure shows that 
the restarted jets are directed towards the brightest warm 
spots at the ends of the relict lobes, indicating that the 
restarted jets are probably tracing the paths of the jets 
just before the central engine turned off. 



2.2 6 cm radio continuum 

The radio image of the inner double source at 6 cm wave- 
length is shown is Figure [2] This image was made with 
the combined BnA and CnB VLA visibilities. The im- 
age has an rms noise of 20 /ijy beam -1 and the dynamic 
range is about 125. This high-resolution 6 cm image does 
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not have the surface brightness sensitivity to detect the 
relatively fainter parts of the inner lobes closer to the 
core, which are traced in the 22 cm image with poorer 
resolution. The emission trail of the NW inner source ex- 
tends approximately 20" from its end towards the core 
and fades out at approximately 55" from the core. The 
emission trail of the SE inner source extends twice as far 
before the surface brightness drops to the same detection 
limit. The 6 cm images have beem made with a circular 
beam of FWHM 1 "3; there are ~ 5 synthesized beams 
across the width of the inner source, and we have well 
resolved all the complex structure within the length of 
the observed inner source as well as at the ends of the 
double structure. 

The inner sources have relatively bright ends and 
lower surface brightness emission trailing off towards the 
core. The emission trails have an almost constant decon- 
volved FWHM about 8"(» 16 kpc). The bright structures 
at the leading ends are narrower and have deconvolved 
FWHM of 3 "-4"{k 6-8 kpc). As observed projected on 
the sky, the bright ends of the inner structures are at 
somewhat different distances from the core: we measure 
the arm-length (distance from the core to the ends) ratio 
to be approximately 1.22. 

Our new high resolution 6 cm wavelength VLA im- 
age reveals a number of characteristics, which are not 
apparent in previous ATCA images of the source. First, 
patchy emission features, which are 2-4 times brighter 
than the underlying emission, are observed along the trail 
in the SE inner source at 6 cm wavelength (see Figure [2}. 
These features are also observed at the same locations in 
a 22 cm wavelength VLA image at 2'.'6 resolution, which 
is shown in Figure indicating that they are genuine. 

The second new feature is that these emission en- 
hancements appear to be located on both sides of a cen- 
tral channel along the length of the trail: there appears 
to be a channel running along the center of the SE inner 
lobe where there is a decrement in the surface bright- 
ness level. A central decrement was also noted in previ- 
ous 13 cm ATCA im ages made with a lower resolution 
jSaripalli et alj|2003l ): however, the relatively poorer res- 
olution had led to the inference that the inner lobe had 
the form of an emission sheath whereas with the improved 
resolution of our new 6 cm image it is apparent that struc- 
ture is that of a cylinder with an emission decrement 
along a narrow cylinder running along the axis. The cen- 
tral decrement in the surface brightness is traced all of the 
way to the head, where it ends at the bright peak. Addi- 
tionally, our new 6 cm image reveals that the decrement 
in surface brightness follows a straight line within the 
southern emission trail and, surprisingly, the decrement 
does not follow the overall curvature that is observed over 
the length of the inner source. A mean profile across the 
width of the SE lobe, which clearly shows the central 
channel, is shown and discussed later in § 14.11 We do 
not observe a central decrement in the surface brightness 
distribution over the short length of the northern trail in 
our VLA 6 cm wavelength image. However, a dip in the 
surface brightness level is detected in profile cuts across 
the faint northern trail in the 22 cm VLA image made 
with 2" 6 FWHM beam (see Figure [3]). The decrements 



trace a straight line from the core to the bright peaks at 
the ends of the inner double. 

The third new feature is that our new VLA 6 cm 
wavelength image resolves the bright ends of the inner 
double radio source. The image reveals extraordinary 
structures at the leading ends, which were not seen in 
previous ATCA images. The SE inner source terminates 
in a rim of enhanced emission, which is approximately 
10 times brighter than the trail towards the core. The 
curved rim extends approximately 5" from the tip along 
the edges of the two sides, and is the brightest feature 
in the radio galaxy, with a peak surface brightness of 0.4 
mjy arcsec -2 . The NW inner lobe terminates in a well re- 
solved, relatively high surface brightness plateau, which is 
~4 times brighter than the trailing emission. The bright 
peaks at the ends of the inner sources are sightly recessed 
from the leading edges of the detected emission, at least 
in projection. 

The total flux density of the inner double source 
at 6 cm wavelength is 35 mjy, implying that the to- 
tal radio power of the inner double source is about 
2.4 x 10 24 W Hz" 1 at 1.4 GHz. The absolute magnitude 
of the host galaxy is M b — —20.9 and, according to the 
relation in lOwenl (| 19931 ). the total radio power of the in- 
ner lobes of the restarting radio galaxy is below the FR 
I-II dividing line. 



2.3 Polarization and rotation measure 

We have made images of the polarized intensity at 6 and 
22 cm wavelengths using, respectively, the VLA CnB 
array visibilities at 6 cm and BnA array visibilities at 
22 cm. The polarization images, at a resolution of 3"5, 
were used to compute distribution of rotation measure 
(RM) over the inner double. We find that the RM is 
fairly uniform, with no significant variations, over each 
of the two components. The NW component has a mean 
RM of —6 rad m~ 2 with a 1-a scatter of 4 rad m~ 2 . 
The SE inner source has a somewhat lower mean RM 
of —15 rad m~ 2 with a 1-a scatter of 5 rad m~ 2 . The 
mean RM over the entire inner double is —11 rad m -2 , 
wh ich is close to t he me an value of —14 rad m -2 found 
bv lSaripalli et all (2003) for the outer lobes. 

In order to examine the polarization structure in 
greater detail over the inner double, polarization images 
with beam FWHM l"3 were made using the BnA and 
CnB visibilities at 6 cm wavelength. In Figure 2] we 
show the distribution of polarized intensity at 6 cm wave- 
length. Overlaid are vectors showing the orientation of 
the projected _E-field, with lengths proportional to the 
fractional polarization. The orientations of the vectors 
have been corrected for Faraday rotation, assuming a 
uniform RM — —11 rad m -2 . Peaks in the polarized 
intensity are observed at the ends of the northern and 
southern sources. 

At the leading end of the SE component, there are 
two peaks in the polarized intensity located to the NE and 
SW of the peak in total intensity. In the trail behind the 
total intensity peak there are two rails of enhanced po- 
larized emission. The fractional polarization is enhanced 
along the boundaries of the source and decreases towards 
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the centre-line where there is also a decrement in the to- 
tal intensity. The patchy emission features that were ob- 
served in total intensity images along the length of the 
SE component are not coincident with peaks in the frac- 
tional polarization. The mean level of fractional polar- 
ization over the SE component is ~ 25%. Around the 
rim in the total intensity image, where there are peaks in 
the polarized intensity, the mean fractional polarization 
is ~ 35%. The projected magnetic field vectors, which 
are perpendicular to the _E-field orientations, are aligned 
with the total intensity contours around the boundaries 
of the source. 

The polarization characteristics of the NW compo- 
nent are similar to those of the SE component. The mean 
fractional polarization over its length is ~ 25%, similar 
to that seen over the SE component. The fractional po- 
larization too is enhanced along the boundaries. There 
are several patches where the fractional polarization is 
relatively low; these are located in regions where the ori- 
entations of the iJ-vectors change sharply in direction, 
and are likely a result of beam depolarization. Along the 
boundaries of the northern source, the magnetic field is 
aligned with the total intensity contours. In the central 
regions of the head and trailing emission, the magnetic 
field vectors are in a direction perpendicular to the radio 
axis. 

We do not find any evidence for Faraday depolariza- 
tion between 22 and 6 cm wavelengths over the NW inner 
source, for which we measure an average depolarization 
ratio (DR; the ratio between the percentage polarization 
at 22 cm to that at 6 cm) of approximately 1. However, 
there is an indication of depolarization over the SE in- 
ner source, where th e average DR is approximately 0.85. 
ISaripalli et ail (|2003h estimated the DR in the outer lobes 
to be in the range 0.96-1.0 at the ends and decreasing to 
approximately 0.92 at the centre, with no DR asymmetry 
between the two lobes. 



2.4 Spectral index 

The distribution of spectral index (S v oc u a ) over the in- 
ner double, between 6 and 22 cm wavelengths, was com- 
puted using VLA visibilities obtained in scaled arrays: 
CnB data at 6 cm and BnA data at 22 cm. Images at 
both wavelengths were made with beams of FWHM 3'.'5. 
Figure [5] shows the resulting spectral index distribution; 
pixels in the individual images that had intensities less 
than 4x the rms image noise were blanked. We observe 
no obvious variations in the distribution of spectral index 
over the two components of the inner double. The NW 
and SE components of the inner double have mean spec- 
tral indices of — 0.65±0.12 and — 0.68±0.15 respectively. 
These values are in agreement with the values found by 
ISaripalli et ail <|2003h between 22 and 12 cm wavelengths; 
therefore, there is no evidence for a curvature in the spec- 
tra of the two components between 22 and 6 cm wave- 
lengths. 

We have examined the data for trends in the mean 
index along the length of the source by averaging the 
spectral index distribution shown in Figure [5] across the 
width of the source. The resulting profile shows no sig- 



nificant trend: we find that along the observable lengths 
of the inner northern and southern sources, from their 
bright ends to their fading emission trails, the spectral 
index remains fairly constant, and at a value in the range 
-0.65 to -0.70. 

We have used our new and improved 22 cm im- 
age of the source, w ith the 13 cm image presented in 
ISaripalli et all (|2003h . to examine the spectral index dis- 
tribution in the outer lobes. We do not include an im- 
age of the computed spectral indices her e, since our re- 
sults are consistent with those derived in ISaripalli et all 
(|2003l ). Overall, there is a steepening of the spectral in- 
dex from the outer extremities towards the centre. The 
spectral indices in the warm spots at the ends of the NW 
and SE outer lobes (a — —0.67 and a = —0.70 respec- 
tively) are very similar to that in the emission peaks in 
the NW and SE inner lobes (a = -0.67). The outer lobe 
regions with higher surface brightness, including those 
regions adjacent to the peaks at the ends of the inner 
double, have relatively flatter spectral indices. The lower 
surface brightness bridge in the vicinity of the core and 
the inner double has a spectral index that is relatively 
steeper and patchy with a in the range —1.0 to —2.0. 



3 THE HOST GALAXY AND 
ENVIRONMENT 

3.1 B and R band images 

We obtained optical images of the host galaxy of 
PKS B1545-321 with the Wide Field Imager (WFI) 
at the prime focus of the 3.9 m Anglo- Australian Tele- 
scope (AAT) in August 2004 during service time. We car- 
ried out observations in the R-, V- and B-bands using 
one of the eight 2k x 4k CCDs. Conditions were non- 
photometric with a mean seeing of 2". The data were 
reduced using standard techniques in the IRAF software 
package. Frames were bias subtracted and flat-fielded us- 
ing dome flats. Image registration was carried out using 
stars in the field, whose coordinates were measured on 
the SuperCOSMOS (SCOS) B-band image with an rms 
accuracy of 0.1 pixels. Reduced images in the R- and B- 
bands are shown in Figures [6] and [7] 

The host galaxy appears to be ordinary in the R- 
band AAT image. Isophotal contours, which are shown 
overlaid on the image, appear to be elliptical and are sym- 
metrical about the centre. In contrast, the host galaxy 
appears somewhat disturbed in the B-band image: there 
is an offset between the centres of the inner and outer 
isophotes and the isophotes are extended to the west. 
Additionally, the B-band image shows a distinct com- 
pact object to the SE and embedded within the faint 
envelope of the host galaxy; however, no corresponding 
component is apparent in the R-band image. It was pre- 
viously noted that the B- band digitized SuperCosmos 
image of the host galaxy showed signatu res of a central 
dust lane perpendicular to the radio axis l)Saripalli et al.l 
2003); however, our new higher quality AAT images show 
no evidence for a dust lane. 

We fitted elliptical isophotes to the R- and B-band 
images using the IRAF routine ELLIPSE (|jedrzeiewskil 
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1987). We ignore the fits to the central 2" diameter re- 
gion since they would be compromised by any ellipticity 
in the point spread function. Immediately outside this 
zone, the fitted elliptical isophotes are similar in the red 
and blue images and have position angles in the range 
60-80°. At a distance of 3.6 kpc (1.8 arcseconds) from 
the nucleus, the fitted isophotes abruptly rotate through 
> —120° and increase in ellipticity. In order to reveal any 
morphological patchiness, we constructed models of the 
host galaxy in the R and B-bands and subtracted these 
from the two AAT images. The IRAF task BMODEL 
was used to create the models by interpolating the ellip- 
tical isophotal fits to the data. We did not include the 
third and fourth order harmonics of the elliptical fits in 
the model subtraction as these were found to introduce 
spurious features in the residual images. 

The residual R- and B-band images are displayed in 
Figures [5] and [DJ Grey-scales in these subtracted images 
signify excess emission over that expected in the models. 
The R-band residual image does not reveal any peculiar 
features within the brighter regions and envelope of the 
host galaxy. On the other hand, the B-band residual im- 
age, clearly reveals a distinct object to the SE of the cen- 
tre and at a projected distance of 10.5 kpc. The counts 
in this region exceed seven times the rms noise in the 
region where the model has been subtracted. Separately, 
there is a hint of another distinct feature to the west of 
the nucleus, at a level of four times the rms noise. These 
components, which are relatively blue in color, may be 
associated with star-forming regions in the host galaxy 
and an indication of a recent merger. 



3.2 The neighborhood of PKS B1545 321 from 
the 6dF Galaxy Survey 

In Figure [TD] we have plotted the positions of all galaxies 
in the vicinity of the host galaxy in redshift space and sky 
posit ion as measured by the 6-degree-Field Galaxy Sur- 
vey (I Jones et al.ll2003 . [20051 ). in a set of redshift slices. 
The 6dF Galaxy Survey (6dFGS) is a near-infrared se- 
lected redshift survey with median redshift z — 0.05, 
complete to K = 12.75 over the entire southern sky 
> 10°). The five panels from (a) to (e) each cover 40°- 
fields (280 Mpc) centered on the host galaxy. Each slice 
is 20 Mpc deep in redshift space, and successive redshift 
slices overlap by 10 Mpc along the line of sight. Panels (a) 
and (b) together cover a depth of 30 Mpc in front of the 
host galaxy; panel (c) covers a depth of ±10 Mpc about 
the host galaxy, and panels (d) and (e) cover a depth 
of 30 Mpc behind the host galaxy. In panels (a) to (e) 
of Figure [TD] we also show in grey scale the local galaxy 
density (per sq. deg.) within the individual redshift inter- 
vals, smoothed by 5°-diameter windows and sampled at 
cells spaced 1° apart. Individual 6dFGS fields are 5?7 in 
diameter and overlap to varying extents. An image show- 
ing the completeness of the survey for the sky region is 
also shown in panel (f) of Figure [TD1 it may be noted here 
that the survey boundary representing the Galactic lat- 
itude limit of the survey lies 7° from the host galaxy of 
PKS B1545— 321. The effects of redshift incompleteness 
have been corrected for in panels (a) to (e) by weight- 



ing the densities using the redshift selection function for 
6dFG S as a function of sky position (see, e.g.. I Jones et al.l 

(Hool)). 

The host galaxy is not in a high over-density. Frac- 
tional over-density factors up to 8 (which is at the 4 sigma 
level) are observed within the 280 x 280 x 41 Mpc cube 
on this smoothing scale. At the location of the host there 
are no 6dF galaxies within R=15 Mpc radius; however, 
the over-density factors within a few degrees on the sky 
and within a few 10's of Mpc take on values up to 3, 
indicating that the host galaxy is not in a void either. 
In velocity space and at redshifts somewhat lower than 
that of the host, there appear to be galaxy over-densities 
to the north and west of the host. At redshifts some- 
what beyond the host galaxy, the over-densities are to 
the south and west. The distribution is complex and an 
elucidation of the detailed structure would require tar- 
geted multi-object spectra of numerous galaxies close to 
the host on the sky, perhaps using the AAOmega on the 
AAT. On the basis of the distribution in 6dF galaxies, we 
believe that the host is located in a relatively low galaxy 
density environment in which the fractional over-density 
is less than 3. It may be noted here that this result is 
on a smoothing scale of R=15 Mpc, and the fractional 
over-density factor could be higher on smaller scales. 

The bottom line is that there is no evidence in the 
6dF distribution that the host is in a rich cluster or rela- 
tively dense parts of the large scale structure of the uni- 
verse; the host is also not in a void. We may conclude that 
the host is located in the relatively low density parts of 
the filamentary galaxy distribution, which has typically 
fractional over-densities in the range 5-200. Assuming 
that galaxies trace the gas in the large scale distribution, 
as is the cas e in large-scale cosmological hydrodynam ical 
simulations (|Cen fe Ostrikerlll999l ; iDave et al.ll200ll ). we 
may expect the IGM gas density in the vicinity of the 
Mpc-scale radio source to be overdense by a factor of 3 
with respect to the mean baryon density in the IGM as- 
sociated with the filaments. Again, here we are assuming 
that there is no small scale structure in the gas in the 
vicinity of PKS B1545-321. Half of the baryons in the 
universe are missing, presum ably in the Warm-Hot IGM 
associated with the filaments l|Cen fe Ostrikeril 19991 ). and 
so the gas density environment of the radio source is at 
least 3 x 0.5 x fWyonPcriticai = 6 x 10~ 31 g cm -3 . 



4 DISCUSSION: THE PHENOMENOLOGY 
ASSOCIATED WITH THE RESTARTING 
JETS IN PKS B1545-321 

4.1 The restarted jets 

On the basis of the Clarke et al. simulations, high resolu- 
tion images of inner doubles in restarting radio galaxies 
are expected to reveal collimated jet structures, with only 
weak hot spots at their ends. Large-sca le jets in power- 
ful radio galaxies are highly collimated (|Bridle fe Perlevl 
1984); typically, t hey have widths o f the order of a few 
kpc {e.g., 3C 353 ISwain et al.lll996l ). However, the ends 
of the inner double structure in PKS B1545— 321 are rel- 
atively bright with deconvolved widths of about 7 kpc. 
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Moreover, the trails of edge-brightened emission between 
the extremities and the core are broader and have de- 
convolved widths of about 13 kpc. It is, therefore, un- 
likely that the entire inner double structure observed in 
PKS B1545— 321 represents collimated restarted jets. 

Within the edge-brightened inner double structure 
we do not observe any relatively narrow, collimated, jet- 
like emission features that may represent emission from 
the restarted jets themselves. We do, however, observe 
a central channel along the length of the southern inner 
component along which there is a decrement in the sur- 
face brightness (see § 12.2)1 . A decrement is also seen along 
the length of the short NS trail, excluding the bright end 
(see Figure [3} . These dips in the surface brightness may 
be the signatures of Doppler dimmed jets. For example, 
jet material moving with a jet Lorentz factor of Fj c t = 7 
on the plane of the sky would have a Doppler dimming 
factor of 0.14. Alternatively, the decrement may be the 
result of the jet being relatively dim because the parti- 
cles in the jet are not reaccelerated and the surrounding 
plasma is freshly accelerated. 

In Figure 1111 we show the mean transverse profile 
across the southern inner source. This slice profile was 
produced by averaging the trailing emission along the 
length of the source, excluding the bright rim of emis- 
sion at the end. The profile shows clearly the central 
decrement in surface brightness. We have modelled this 
profile as representing emission from an optically thin, 
thick- walled tube, in which the material in the tube wall 
has a constant emissivity. We convolved the 1-D slice ex- 
pected from this model with the synthesized beam and 
then fitted the resulting function to the observed profile 
data. The best fit (also shown in Fig. 1 1 1 j) corresponds 
to an outer diameter of 20 kpc for the cylinder and a 
diameter of 3.2 kpc for the inner hollow section. Our hy- 
pothesis is that the central hollow in the inner lobes of 
PKS B1545— 321 represents Doppler-dimmed restarting 
jets. The observations indicate that the inner lobes are 
not jets, but composed of central jets surrounded by in- 
ner lobes, which in turn are embedded within the relict 
outer lobes. 

These jets on the two sides of the core appear to be 
misaligned (i.e. the inferred jet paths on each side of the 
core are not collinear) . Moreover, the jets appear to trace 
a straight path to the ends of the inner double whereas 
the inner lobes show mild curvature with a C-shape or 
reflection symmetry. The structure is not inversion sym- 
metric, as we might expect if the jets are precessing. Pre- 
sumably, the restarted jets were initially collinear and 
have gradually evolved to their present misaligned state. 



4.2 Are there bow shocks in the relict cocoon? 

Bow shocks are exp ected to lead the advance of restarte d 
jets in old cocoons (|Clarke fc Burnslll99ll ; IClarkelll997n . 
We expect to observe these shocks via enhanced syn- 
chrotron emission since the relict lobes of a radio galaxy 
contain relativistic particles and magnetic fields, and the 
bow shocks should enhance the emissivity by compres- 
sion of the relativistic plasma and reacceleration of the 
electrons. In the synthetic surface brightness images de- 



rived from the restarting jet simulations, the leading bow 
shocks do produce observable feature s that are d istinct 
from the relict outer cocoon emission (|Clarkell 19971) . The 
bow shock would be expected to weaken with distance 
from the apex, but might remain an observable feature 
over its entire length. At the apex, where the bow shock is 
stronger and balances the terminal shock of the new jet, 
the brightness contrast might be an order in magnitude 
or more. In sharp contrast to the predictions of the nu- 
merical simulations, bow shocks ahead of restarting jets 
have not been evident i n observations of double-d ouble 
radio sources to date fe.g. lSaripalli et al! 2002, 2003]), and 
this has been a major motivation for the high resolution 
and sensitivity radio observations presented here. 

In our 22 cm image (Fig. [TJ we measure a mean sur- 
face brightness of 450 /iJy beam -1 and 600 /iJy beam -1 
in the relict cocoon emission in the vicinity of the ends 
of the NW and SE inner double; the image rms noise is 
40 fj,3y beam -1 . Within the errors in the image, there is 
no evidence for wide bow shock like features in the outer 
relict cocoon exterior to the inner double. While there 
exist regions of enhanced surface brightness around the 
inner double lobes, notably along the western side of the 
NW inner component, these regions are patchy and may 
be extensions of the relatively brighter material that is 
observed ahead of the inner double lobes. If there is a 
wider bow shock, which leads the advance of the new 
jets in the relict lobes, then we estimate that the bright- 
ness of any such bow shock, relative to the cocoon is 
less than about 20%. Relative to the hot spots at the 
ends of the restarted jet, the brightness of such a bow 
shock is less than 10%. Our improved images confirm 
that PKS B1545— 321 does not have bow shocks exterior 
to the inner double and with a br ightness expected from 
simulations (in, e.g., lClarke| ll997l '). 

However, we note here that the bow shock may be 
weak over most of its length and, therefore, reacceleration 
of particles in the relict cocoon plasma may be insignif- 
icant. In this case, any brightness enhancement as a re- 
sult of the bow shock would be the result of compression 
of the relict lobe plasma. If the relict cocoon emissivity 
has a straight synchrotron spectrum with spectral index 
a = — 1, and a tangled magnetic field, 1-D compression 
by factor / by a bow shock would enhance on the emissiv- 
ity by factor / 8/ ^ 3 : the limit of 20% on the enhancement 
of emissivity, which we have derived above, requires that 
the compression be at most 7%. Curvature in the cocoon 
electron energy spectrum, which might be expected in the 
relict plasma as a result of enhanced radiative cooling of 
the more energetic particles, may weaken the constraint 
on the bow shock compression ratio. However, the ex- 
tremely low compression ratio required by the stringent 
observational limits placed by the data presented herein 
is unlikely. 

4.3 The hot spots at the ends of the inner lobes 

4.. 3.1 Total intensity features 

Bright hot spots at the ends of radio lobes represent 
the working surface where the jet meets the ambient 
medium. Usually this medium is the relatively dense 
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inter-stellar medium (ISM) or IGM, and the result is 
a strong shock, which converts a large fraction of the 
jet kinetic energy into relativistic particles. The ambient 
medium of the restarted jets of PKS B1545— 321 are the 
relic synchrotron lobes of the giant radio source. These 
are expected to be lighter than the ISM or IGM, even if 
there has been entrainment of thermal material into the 
lobes. For this reason, we might expect the appearance 
of hotspots at the ends of restarted jets to be unusual. 
In addition, we might also expect that the hotspots at 
the ends of restarted jets are a merger of both the termi- 
nal jet shocks and the bow shocks immediately ahead of 
the jets since both shocks should reaccelerate relativistic 
particles and amplify magnetic field. 

The bright ends of the inner double source of 
PKS B 1545— 321 are resolved in the 6 cm VLA image 
(Fig. EJ) . The SE source terminates in a bright extended 
peak, where as the NW source terminates in a broad 
plateau of emission. The peak at the end of the SE source 
resembles hot spots created by powerful jets in the IGM; 
however, the end of the NW i nner double does no t show 
such a bright compact feature. iLeahv et alJ (| 19971 } define 
a hot spot as any feature that is not part of a jet, that 
has a largest dimension smaller than 10% of the main 
axis of the source, a peak brightness that is greater then 
10 times the rms noise in the image, and that is separated 
from neighboring peaks by a minimum, which drops to 
less than two-thirds of the brightness of the fainter peak. 
If we adopt this definition, then the peak at the tip of the 
SE inner source may be classified as a hot spot. By fit- 
ting a Gaussian to the peak at the end of the jet, we have 
estimated the intrinsic size of the major and minor axes 
of the central compact feature to be 3'.'l and 2 '.'1 respec- 
tively, corresponding to linear sizes of 6.0 and 4.1 kpc. 
The measur ed size is consistent with the size distribution 
observed bv lHardcastle et alJ |l998) for a sample of pow- 
erful z < 0.3 radio galaxies. As far as the NW structure is 
concerned it would be overprescriptive to insist that this 
region strictly satisfy the Leahy et al. hot spot criteria. 
The existence of a region at the NW end of this structure, 
which is significantly brighter than the rest of the tail, is 
at least consistent with the interaction of a new jet with 
the pre-existing lobe. Also, in view of our interpretation 
of the arm-length asymmetry of the inner double being 
due to relativistic motion, this region would be Doppler 
dimmed (see § I5.lt . 

Despite the similarity to classical hot spots, we note 
that the compact feature at the end of the SE inner lobe 
displays some unusual characteristics that are atypical of 
hot spots observed at the extremities of powerful radio 
galaxies and are, on the other hand, what we might ex- 
pect of jets developing in a low density environment. In 
particular this feature shows an extended rim of emission, 
which wraps around the head of the source. This rim is 
extended in a direction perpendicular to the radio axis 
and shows tails, which extend on either side of the cen- 
tral jet (identified to be the central decrement). The tails 
extend approximately 5" (10 kpc) back toward the core of 
the radio galaxy. We have examined high resolution im- 
ages of powerful radio galaxies in the literature, and do 
not find any hot spot that resembles the compact struc- 
ture at the end of the SE inner lobe; no hot spot observed 



to date shows such a neat rim wrapping around an incom- 
ing jet. Indeed many hot spots project beyo n d the envi- 
sion of the associated lobe (see iPerlev et al.1 l|l997l ) and 
references therein). The ends of the lobes of powerful ra- 
dio galaxies often show complex structures, and in many 
cases multiple hot spots. Although there exist examples 
of multiple hot spots distributed in rim-like structures at 
the ends of sources (e .g., the northern lobe of 3C173.1, 
lHardcastle et al1ll997T ) it is difficult to find individual hot 
spots with swept-back wings or tails. The western lobe of 
3C234 has a central compact compone nt and ridges on 
each side that bend back into the lobe |Hardcastle et all 
1997); however, in this case the northern ridge is likely 
to be a distinct and separate hot spot. We suggest that 
the extended rim-like structures represents a blending of 
emission from the jet terminal shocks and the immedi- 
ately preceding bow shock. 

The compact feature at the end of the SE inner dou- 
ble is the brightest feature in PKS B1545— 321. However, 
it has a peak surface brightness of only 0.4 mjy arcsec -2 
at 6 cm wavelength, which is significantly smaller than 
the typical surface brightness observed for hot spots 
in powerful radio galaxies. For examp le, in the sample 
of pow erful radio galaxies presented in lHardcastle et al.l 
(1 19971 ) the surface brightness values of the hot spots are 
an order of magnitude or more higher. The terminal fea- 
tures at the ends of both the inner SE and NW lobes 
are intermediate in brightness between the strong hot 
spots observed in FR2 radio sources and the low con- 
trast Mach disk structures formed at the ends of the 
new jets in the hydrodynamic resta rting jet simulations 
l|Clarke fc Burnslll99ll ; lciarkelll997T ). We further suggest 
that this is indicative of jets terminating in an environ- 
ment intermediate in density between the normal IGM 
and the extremely low density of a lobe which has not en- 
trained significant external gas. This assertion is treated 
quantitatively below. 



4-3.2 Polarization features 

The polarization characteristics of the inner components 
also provide useful information. At the SE termination, 
the magnetic field runs around the head of the source, 
neatly following the tails in the total intensity distribu- 
tion (see Fig. [4}. At the NW end, the magnetic field is 
aligned with the total intensity contours. In the central 
regions of the broad plateau the magnetic field is trans- 
verse to the source axis. The observed polarization prop- 
erties at the ends of the inner northern and southern 
sources are consistent with what is observed in the hot 
spots of powerful radio galaxies, where the magnetic field 
is nearly always perpendicular to the source axis. Using 
images with 3'.'5 FWHM beam, we measure the spectral 
index at the ends of the NW and SE inner sources to be 
a ~ —0.7 (see Fig. [5}. This is also consistent with that 
observed in the hot spots of powerful radio galaxies. 

Notwithstanding the similarities with the polariza- 
tion properties of classical hot spots there are some in- 
triguing differences. For exam ple the hot spots i n classic 
FR2 sources such as Pictor A (|Perlev et al.lll997n exhibit 
strong polarization coincident with the brightest part of 
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the hot spot with the position angle of the E- vector indi- 
cating a magnetic field perpendicular to the jet direction. 
In PKS B1545— 321 there are local minima in the po- 
larized intensity either coincident with or nearby to the 
brightness maxima. The directions of the E-vectors along 
the boundaries of the SE and NW inner components in 
PKS B1545— 321 are consistent with the expected align- 
ment of the magnetic field along a bow shock. The min- 
ima near the positions of peak intensity are consistent 
with the superposition of comparably bright regions in 
which the magnetic field is (a) perpendicular to the jet 
(the main hot spot emission) and (b) aligned almost par- 
allel to the jet (the bow shock emission). There are a 
number of local maxima and minima in the polariza- 
tion images which could possibly be interpreted as the 
result of jittering of the jet near the respective hot spots. 
However, in view of the straightness of the jets inferred 
from the transverese intensity profiles, we feel that not 
a great deal of jittering is occurring in this region of the 
source. On the other hand the local maxima could indi- 
cate regions where the inner bow-shocked material inter- 
acts with inhomogeneities in the cocoon. 



5 ESTIMATES OF SOURCE PARAMETERS 

We have noted that there are several possible restart- 
ing jet models, which may account for the inner struc- 
tures observed in PKS B1545— 321. In the sections be- 
low, we derive estimates for a number of physical param- 
eters associated with the source, including the density of 
entrained matter in the relict cocoon, and the external 
Mach number of the ends of the restarted jets. These es- 
timates allow us to derive some conclusions about the 
formation mechanism of the inner double in this source. 



5.1 Arm length and brightness asymmetry in 
the inner double 

Consideration of the arm length asymmetry in this ra- 
dio galaxy provides a valuable constraint on the physical 
ideas discussed in section [T] One expects the lobes of 
a powerful radio galaxy to be significantly lighter than 
the ambient ISM or IGM. Hydrodynamical simulations 
of restarting jets show that in case of an interruption, 
restarted jets are of comparable density to the cocoons 
of the old jets, and as a result they propaga te rapidly 
through the old cocoons l|Clarke fc Burn j| 199 ll ). As com- 
pared to the original jets in PKS B1545— 321, we expect 
that the new jets would propagate relatively unimpeded 
in their advance from the core, unless there has been sig- 
nificant entrainment of thermal material into the relict 
lobes. If the twin oppositely directed jets are advancing 
at equal and constant velocities, then time retardation 
introduces an asymmetry in their apparent projected dis- 
tances from the core; an argument for constant advance 
velocities is given in section HJ] In the usual case where the 
jets terminate in the ISM or IGM, caution is required in 
interpreting arm-length ratios since asymmetric distribu- 
tions of matter around the galaxy can overwhelm time re- 
tardation effects. However, in PKS B1545— 321 the local 



ISM/IGM has been cleared and the new jets are propa- 
gating into synchrotron bubbles created by past activity. 
This is probably the most ideal circumstance under which 
time retardation effects may be apparent. This would be 
the case even if there is entrainment of the local IGM, 
since the entrainment rates should be similar on both 
sides of the source. If there is any unevenness in the en- 
trainment rate (as, for example, a result of inhomogeneity 
in the local IGM distribution) , this should be evened out 
by the dispersal of entrained gas into the large volume of 
the cocoon. 

Let dhead and d Ci hcad be the projected lengths of the 
ends of the jet and counter jet as measured from the core, 
/3hs the speed of the hot spots relative to the speed of light 
and 9 the angle of inclination of the jet axis to the line 
of sight. Then the arm length ratio is 

dhoad 1 + Ais cos 9 



^c.hcad 1 

implying that 

1 D- 1 



f3 hs cos ( 



Phs = 



cos e D + 1 ' 



(1) 



(2) 



The observed arm length ratio for the inner double in 
PKS B1545-321 is D = 1.22 ± 0.02. The corresponding 
values of /3hs are plotted in Figure[l2]for angles of inclina- 
tion in the range 45° < 9 < 90°. Since PKS B1545-321 is 
a radio galaxy, and not a quasar, we expect that 9 > 45° . 
An immediate implication of equation [2] is that the the 
angle of inclination cannot exceed 9 = 84? 3, since the 
hot spot speed /3hs < 1. The plot shows that for the 
allowed range of possible viewing angles (45° — 84? 3), 
/3hs takes on values in the range 0.14-1.0. The lower 
limit of about 0.14c on the advance speed of the new 
jet is an order of magnitude larger than the typical ad- 
vance spee ds inferred for powerful jets ploughing into the 
ISM/IGM (jScheuerl 19951 ). consistent with the model that 
the restarted jets in PKS B1545— 321 are indeed evolving 
into a lighter than usual medium. 

The brightness asymmetry ratio in hot spots is an- 
other consequence of relativistic speeds at the ends of 
jets; the observed ratio of the brightness in the hot spots 
at the NW and SE ends of the inner double is approx- 
imately 2. The brightness asymmetry is related to the 
arm-length asymmetry by 



R = D 



(3) 



where a is the spectral index of the hot spots. Given our 
observed values of a = —0.7 and D = 1.22, the brightness 
asymmetry is expected to be R ~ 2.1, which is appeal- 
ingly consistent with our measured value of 2. 

The measured arm- length asymmetry ratio D = 1.22 
implies that /3h s cos# ~ 0.1. Giant radio sources are 
selected from compilations of large angular size radio 
sources. Hence, are biased towards galaxies lying close to 
the sky plane with relatively larger angles of inclination. 
This suggests larger values for /3h s • At extreme angles of 
inclination where /3h s takes values close to 1, we expect 
Ais ~ /3jet , and a ballistic advance of the new jets through 
the outer synchrotron lobes. Thi s is the case predicted 
by the resta rting jet simulations l|Clarke fc Burns] Il99ll ; 
IClarkdll997h . However, the range of angles for which /3hs 
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approaches 1 is small. For a large range of viewing angles, 
it is likely that /3h s < Pjet , in which case the advance is not 
ballistic. This latter scenario is consistent with the mor- 
phology of the inner lobes and relatively bright hotspots 
at their ends. 

5.2 Density of entrained matter in the relict 
cocoon 

Assuming a non-ballistic advance for the jet ends, we 
may estimate the density of the outer lobes by relating 
the hot spot pressure to the cocoon density and the hot 
spot advance speed. In this case, the density of the cocoon 
material, p c , is given by, 

— ^ Phs ,.n 

Pc - 4(r hs /3 ha c) 2 ' w 

where p ns is the pressure in the hot spot and IYs = 
— /3jJ s is the corresponding Lorentz factor. (This ex- 
pression does not account for all relativistic effects, which 
we consider below.) 

We estimate a minimum hot spot pressure of p m i n = 
1.5 x 10 -11 dyne cm -2 at the end of t he SE inner dou- 
ble, using minimum energy conditions l|Milevlll980l ). We 
expect the true hot spot pressure phs -2 Pmin- Allowing 
for relativistic effects, the pressure scales as 5~ 2 , where 
5 — (1 — /3h s cos^)" 1 is the Doppler factor. Therefore 
allowing for relativistic effects in the minimum pressure 
estimate, the density in the relict cocoon ahead of the 
restarted jet is: 



_ 3 Phs Pmin (1 - Phs COS 6) 2 



4p n 



PL 



(5) 



The density implied by equation ([5]) is plotted in 
Figure [13] as a function of angle of inclination 6 for 
Phs/Pmin=l. The relationship between /3 ns and 0, as im- 
plied by the arm-length asymmetry, has been used. The 
plot shows that the cocoon density is in the narrow 
range (0.1 — 5) x 10~ 31 g cm -3 for the range of plausi- 
ble viewing angles. Repeating the calculation for the end 
of the NW trail, which has a minimum energy pressure 
Pmin — 1.0 x 10 -11 dyne cm 2 , gives a similar range of 
density values. 

As shown in § 13.21 we expect that pigm is greater 
than or of the order 6 x 10~ 31 g cm -3 . Hence, the cocoon 
density as a fraction of the IGM density k = p c / pigm is in 
the range of 0.02 - 0.8. Smaller values of re are implied for 
larger angles of inclination. Our derived values of « are, 
once again, consistent with the notion that the restarted 
jets in PKS B1545— 321 are evolving into a medium that 
is lighter than the IGM. At large angles of inclination, 
the cocoon density is small (less than 2% of the IGM 
density). It is under these conditions that we may ex- 
pect the density of the restarted jets relative to the den- 
sity of the external cocoon (77 = Pj e t/p c ) to be of order 
unity or greater, as predicted by the resta rted jet simula- 
tions l|Clarke fc Burnslll99ll ; IClarkelll997h . Alternatively, 
at smaller angles of inclination, where the cocoon den- 
sities represent larger fractions of the IGM densities, we 
expec t rj < 1 and the entrainment model of iKaiser et al.l 
(2000) to be more relevant. We return to this point again 
in SlO 



5.3 Expansion timescales for the outer hot 
spots 

The jets in the powerful radio galaxy PKS B1545— 321 
no longer supply energy to the outer lobes of the giant 
source; they have stopped, restarted, and the new jets 
currently terminate at locations that are well recessed 
from the ends of the outer lobes. As discussed in sec- 
tion 12.11 the presence of warm spots at the ends of the 
outer lobes suggests that the sound crossing times in the 
lobes are insifficient to completely dissipate the evidence 
of the previous jet termination. The time required for hot 
spots with an initial radius Ro and pressure po to expand 
to the observed pressure p ns is given by 



t 



Phs 

po 



-3/4 



1 



PO 



1/2 



(6) 



(see Appendix [Bj . 

Assumin g standard minimum energy conditions 
l|Milevl Il980l ). we estimate the pressure in the remnant 
hot spots, which have a diameter of ~ 15 kpc, to be 



Phs 



x 10 _ 13 dyne cm . Based on the em pirical rela- 



tionship derived by lHardcastle et all i|l998ft . we assume 
that the hot spots at the ends of the 1.05 Mpc pow- 
erful giant radio galaxy had a radius of Ro = 2.5 kpc 
while the jets were actively feeding them, implying (us- 
ing equation IB3|) that the active hot spot pressure was 
po = 8 x 10 -11 dyne cm -2 . Considering the external co- 
coon densities derived in section 15.21 it follows that the 
jets ceased feeding the hot spots at most 1.8 x 10 5 yr ago. 



5.4 The Mach number of the advance of the 
inner double 

5.4-.1 Dynamical estimate 

We may relate the Mach number of advance of the ends 
of the inner double to the speed of the hot spots and the 
parameters of the synchrotron plasma in the relict cocoon 
traversed by the restarted jets. We use the following no- 
tation: p SlC — c SiC /c, where c s , c is the sound speed in 
the relic coco on; r s , c is the corresponding Lorentz factor 
l|Konigjll980h . The Mach number of the new hot spots 
with respect to the sound speed in the relict cocoon is 
given by 



A^hs — rhs/3hs/r s ,c/?s,c- 



(7) 



We use the parameter \ to denote the ratio 
which parametrizes the rest mass density in a medium 
whe re the pressure is dominated by relativistic particles 
(see lBicknellll994h . In the relict cocoon, the density is 
dominated by thermal matter, and the pressure, p c , is 
dominated by the pressure of the relativistic gas, so that 
Xc = p c c 2 /4p c . Similarly, in the jet, xjet = Pjctc 2 /4pj e t. 
From equation ([5} one has 

3 Phs Pmin (1 — /3 hs COs6>) 2 



16 Pmin Pc 



PL 



(8) 



In terms of \c, the sound speed and corresponding 4- 
velocity in the cocoon are given by 



p; 



s^i + xc)- 1 



(9) 
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r s , c /3 s , c = (2 + 3 X c)" 1/2 (10) 

Combining equations ((2} for the velocity in terms of 
the arm- length asymmetry, §8§ for the density parameter, 
and equation J7J for the hot spot Mach number, enables 
us to determine the latter as a function of the inclination 
angle. 

We can use the above relationships to gain a semi- 
quantitative idea of the implications of the inferred pa- 
rameters when Ths ~ 1 and r s , c « 1. In this non- 
relativistic approximation, the velocity dependence can- 
cels out and the Mach number of the new hot spot is 
simply given by: 

M hs « \ (^) (11) 

Assumin g standard minimum energy conditions 
(|Milevl Il980h . the relict cocoon in the vicinity of the 
ends of the inner double has a pressure p c = 3 x 
10 -13 dyne cm -2 . We may expect the true pressure in 
the radio lobe to be close to that given by the mini- 
mum energy condition, since the run of p c versus the 
ratio of particle to magnetic pressures (p par t/(-B 2 /87r)) 
has a shallow minimum, and there is observational evi- 
dence that the conditions in many radio lobes and even 
hot spots are clo se to equipartition l|Croston et alj|2005l : 
lHardcastldl2003l ). Hence, using this value for the cocoon 
pressure and the previously derived value for the hot spot 
pressure, we derive an estimate of Mhs ~ 5 for the hot 
spot Mach number, independent of angle of inclination. 

Let us now consider the more general case where the 
ends of the new jets could advance with relativistic veloc- 
ities. As we have seen in sect ion I5TT1 this is especially rele- 
vant if the radio axis of the giant source is closely aligned 
with the plane of the sky. In this case, the external Mach 
number of the hot spot is unbounded as /3hs approaches 1. 
In Figure [T4l we show the external hot spot Mach num- 
ber as a function of viewing angle, as inferred from the 
above equations for Phs/Pmin = 1 and 10. These curves 
demonstrate that for Phs/Pmin = 1, hypersonic external 
Mach numbers are only implied for a very narrow range 
of viewing angles. For most viewing angles, the external 
hot spot Mach number Mhs ~ 5, and is only supersonic 
rather than hypersonic. On the other hand, the external 
Mach number is higher if the pressure in the hot spot 
is greater than the minimum energy estimate. For exam- 
ple, the phs = lOpmin curve shows that over most angles 
Mhs ~ 15. Hence the inference of a merely supersonic or 
hypersonic Mach number depends upon the ratio of the 
hot spot pressure to the minimum energy pressure. 

5.4-2 Mach angle estimate 

In the hypersonic restarting jet model, the outer edges of 
the inner lobes in PKS B1545— 321 represent the 'missing' 
bow shocks. Thus, it is also of interest to derive an esti- 
mate of the external Mach number of the new hot spots, 
which is based on the opening angles of the putative bow 
shock. 

If the axis of the inner jet is unchanging and the 
inner lobe structures represent relict lobe material that 
has been rejuvenated by bow shocks, then the widths of 



the emission trails are representative of the external Mach 
numbers of the restarted jets. The asymptotic angle of 
the bow shock with respect to the jet axis, #bs, is related 
to the external Mach number of the hot spot by #bs = 
sin" 1 AC 

The trail of the inner Southern double is long enough 
that we can make a reasonable estimate of its inclina- 
tion angle with respect to the radio axis. In order to 
measure the increase in width with distance along the 
southern inner trail, we first rotated the 6 cm image 
by — 33?2 to make the jet axis vertical and then con- 
structed profiles at different distances from the ends, av- 
eraging over axial distances of 2.5 synthesized beams in 
each case. The FWHM of the profiles were measured 
by including all intensity values exceeding 3 times the 
rms noise. The measured half-widths along the SE inner 
source are plotted in Fig 1151 The expression for the bow 
shock angle is an asymptotic one. Therefore we did not 
construct FWHM profiles at the bright ends where the 
width changes abruptly with length. We did not carry 
out the same analysis on the NE trail, since it is much 
shorter and fainter. In order to measure the opening an- 
gles of the low-surface-brightness SE emission trail, we 
least-squares fitted a straight line to the half- width versus 
distance data: The slope of the fit to the SE inner lobe 
data implies a half-opening angle of 2? 5. Simply inter- 
preted, the measured opening angle would imply that the 
SE hot spot is advancing with a Mach number Mhs ~ 23 
with respect to the gas in the old cocoon. Thus there is 
a discrepancy between this Mach angle-based estimate of 
Mach number and the dynamically based estimate ~ 5 
derived above in § 15.4.11 

There are at least three possible explanations for 
this. One possibility is that the jet changes direction 
slightly over its lifetime so that the end of the lobe is 
wider than that associated with a jet which remained 
constant in direction. This would decrease the apparent 
shock angle. Another possible explanation is that the in- 
ner lobes are affected by inhomogeneities in the cocoon 
resulting from the entrainment of IGM. In this case, one 
expects the density to increase towards the sides of the 
coccon and this would have the effect of keeping the 
inner lobes narrow. A third possibility is that the hot 
spot pressure is more than an order of magnitude greater 
than the minimum ene rgy pressure. Given the results of 
lHardcastle et all l|2004 ) who showed, in a sample of 65 
hot spots, that departures form equipartition are unlikely, 
this possibility is remote, with the proviso that hot spot 
dynamics within a pre-existing cocoon could be different. 

5.5 The energy budget 

The energy budgets corresponding to the original and 
restarted jets provide additional information on the ages 
and advance speeds of the outer and inner radio lobes. 

5.5.1 The inner lobes 

The FWHM w 3'.'l of the SE hot spot is larger than the 
jet diameter « l'.'6, inferred from the surface brightness 
decrement in the SE lobe. This may be the result of the 
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dynamics of the shocked plasma near the jet terminus or 
some wandering about of the jet termin us as envisaged 
in the Dentist Drill model (|Scheuerll 19821 ). Whatever the 
reason, we take the area ^4hs of the hot spot as defining 
the area of the working surface. The balance between the 
momentum in the jet and that transfered to the relict 
cocoon over Ahs leads to the following expression for the 
energy flux in the jet (derived in Appendix [Cj: 



1+Xjo 



xp c c 3 /3 je tr hs 2 



[1 - (/Wftct)] 2 



Ah.. (12) 



The hot spot area A^s ~ tt(3 kpc) 2 . In Figure [TBI we 
have plotted Fe as a function of inclination angle, using 
the arm- length asymmetry relation (equation [2} and the 
external (cocoon) densities (derived from equation(5]l . We 
have also assumed that /3j ct w 1. This plot shows that for 
the plausible range of viewing angles, Fe > 10 44 erg s _1 ; 
for a restricted range of inclination angles « 3.5° in which 
the jet is close to the plane of the sky Fe might be 1-2 
orders of magnitude larger than 10 44 erg s _1 . However, 
this is clearly unlikely. 

We can check the above estimate of the energy flux 
against the energy deposited in the inner lobes. The total 
energy deposited into the inner lobes by the restarted jets 
over time At is, 



U « f ad F E At, 



(13) 



where the factor / at j allows for adiabatic losses; we adopt 
a nominal value of / a( j = 0.5, which is typical of a 
non-thermal bubbl e inflated by a relativistic jet (e.g. 
iBicknell et ail 1 1997ft . 

The SE inner lobe, which has a projected length of 
~ 70 kpc in our VLA 6 cm image, has a minimum energy 
of ?7min = 2.8 x 10 57 erg. This leads to another expression 
for the jet energy flux: 

Fe/10 44 erg s" 1 = 8(3 hs sin6 (14) 

where we have used the relation At = Z b s /AisCsin# with 
^obs = 70 kpc. This equation is also plotted as a function 
of inclination angle in Figure [161 The agreement between 
the two independent estimates of the energy flux is excel- 
lent for inclination angles 6 < 80°. This is independent 
confirmation of our general assumption that the source 
is not within 10° of the plane of the sky. 

5.5.2 The outer lobes 

The total minimum energy in the outer lobes of the Mpc- 
scale radio galaxy is U ~ 6 x 10 59 erg. If we assume that 
the energy flux of the original jets is the same as that for 
the restarted jets, then we can also estimate the active 
lifetime and speed of advance of the old jets. Adopting 
/ad = 1/2 and assuming that the original jets also have 
an energy flux of about _Fe ~ 1-6 x 10 44 erg s _1 it follows 
that the old jets were active for 0. 3-2 x 10 s yr. Assuming 
that the lobes have not expanded significantly in the time 
since the first jets were switched off, this implies an ad- 
vance speed in the range /3hs ~ 0.01-0.06/ sin (9. 

The advance speed of the original jets is therefore 
similar to the ad vance velocitie s inferred for jets in power- 
ful radio sources l|Scheuerl 19951 ): the similarity in advance 



speeds for giants and smaller size sources suggest that the 
effective hot spot areas are larger in giant sources, where 
the hot spots are farther away from the AGNs, compen- 
sating for the lower density media encountered by the jets 
in these sources. 

Finally, we note that if the energy fluxes and effi- 
ciencies of the old and restarted jets are similar, then our 
analysis of the energy budget is consistent with an or- 
der of magnitude difference in hot spot speed in the two 
phases of radio galaxy evolution. 

5.6 Restarting timescales 

Suppose the first jet in PKS B1545— 321 turns off at t — 
ti; then the hot spot at the extremity of the outer lobe 
is fed for a further time L/c/3j ot , where L is the length 
of the old lobe. The hot spot then expands for a further 
time t c 



t = ti + 



Hence the current time is 



+ t c 



cf3 ic 



(15) 



Now suppose that the new jet turns on at t — ti and 
that the hot spot propagates for a distance I at speed 
c/3 hs . This gives a second expression for the current time: 

(16) 



t = t 2 + 



C/Shs 



Hence, equating these two expressions 
_L l_ 

C/?hs 



At = t 2 - h = 



CP- 



+ t c 



'jet 



(17) 



The condition that ti > ti leads to the following condi- 
tion on (5 hs , 

0b. > - -j (18) 

Ctexp + -k/Pjet 

In § 15.31 we computed t cxp < 1.8 x 10 5 yr for the 
relict hot spots, based on the external lobe densities de- 
rived in § 15.21 This value, together with I « 150 kpc, 
L ~ 490 kpc (assuming orientation near the plane of the 
sky) and /3j ct = 1, gives /3h s > 0.27. This suggests that 
PKS B1545— 321 is observed at an inclination angle of at 
least 70° to the line of sight. In light of this result, to- 
gether with the unreasonably high energy fluxes implied 
for extreme inclination angles (see § 15. 5[) . we assume that 
PKS B1545-321 is at an inclination angle 70° < 9 < 80°, 
for the remainder of the discussion in this paper. 

A final point on timescales relates to the time be- 
tween cessation and renewal of jet activity, At. Inserting 
numbers in equation (|17[l for At and taking /3j c t — 1 we 
have 



At = 3.2 x 10 J 



490 



150 



+ texp yr 



(19) 



From our estimate of t cxp < 1.8 x 10° yr, it follows that 
At « (5 — 10) x 10 5 yr. Considering our estimate for the 
age of the Mpc-scale source in § 15.5.21 this implies that 
At is at most a few percent of the source age. 



5.7 The jet Lorentz factor 

In this section we derive estimates of the jet Lorentz fac- 
tor using the above estimates of hot spot velocity. First, 
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(21) 



we assume that Xjot ~ 1 as inferred in a number of ra- 
dio jet sources {e.g. M87. IBicknell fc Beeelmanl Il996l . 
and Markarian 501, IBicknell et al.l 2001). The value of 
this parameter in the cocoon, \ c takes values between 
20-90. This range is computed from the range of co- 
coon densities derived in section 15.21 for 70° < 9 < 80° , 
along with the pressure at minimum energy conditions, 
p c = 3 x 10" 13 dyn cm 2 , in the external lobe. 

For a relativistic jet propagating in a cocoon whose 
pressure is dominated by relativistic plasma, the advance 
velocity of the hot spots is given by 

l/2p a 

As ~ i + ^r jct (20) 

where 

CL+Xjet) / Pjct \ _ / pjetc 2 + 4p jct 
(1 + Xc) V Pc J V PcC 2 + 4p c 

This e quation is equivalent to equation (5) of lMarti et al.l 
(1997). Equation (|20[) is relevant when the jet is straight 
and its momentum flux is not spread over a larger area 
than the jet cross-section. This expression is therefore 
pertinent to the hypersonic jet model. 

In order to use equation (I20|l to determine the jet 
Lorentz factor, we solve for F to find: 

Tlct «) (22) 

In equation (|21l) for the jet density ratio, we cannot di- 
rectly measure the pressure in the restarted jets. How- 
ever, we assume that the restarted jets are in pressure 
equilibrium with the inner lobes, which have a minimum 
energy pressure of 5 x 10 -12 dyn cm -2 . Thus, using the 
above estimates, we derive values of r\ in the range 0.3-2; 
this range is consistent with the notion that the jets are 
evolving into a medium of similar density. 

The inferred range of r\ values, together with equa- 
tion (|22l) for the jet Lorentz factor, imply that Pjet = 2— 
3. Specifically, this if the restarted jets advance with hot 
spot velocities in the range 0.3 < /3hs < 0.6, which are 
relevant to the likely range of inclination angles. 

5.8 Entrainment 

Here we consider the potential origin of the gas inside the 
relict cocoon and show that most of the gas must have 
been entrained. First we estimate the mass density of the 
non-thermal plasma, p nt , that is deposited by the original 
jet into the outer lobes. This depends on the mass flux of 
the original jet, M, the time for which the first jets were 
active, i on , and the volume of the pre-existing cocoon V c : 



Put 



Mtp 



(23) 



In this expression M = PjctPjct^jct^ljct- 

We assume that the original jets have the same 
mass flux as the restarted jets. Using our inferred in- 
ner jet diameter dj ct = 3.2 kpc (see bj4.1[l and an inner 
jet pressure that is equal to the minimum energy inner 
lobe pressure, pj ct = 5 x 10 -12 dyn cm -2 , we estimate 
M = 5 X 10 22 Xjct(r 2 - 1) 1/2 g s _1 for the restarted jets, 
and therefore for the original jets. 



Considering the one-sided volume of the outer ra- 
dio lobes and an active timescale of i on = 1 x 10 8 yr 
(see 3S3Ji we obtain the following expression for the 
non-thermal mass density in the outer radio lobes, 



te«4x io- 34 Xj et(r 2 e 



i\l/2 -3 

1) ' g cm 



Adopting Xje 



1 and r jet 



2 gives p n 



(24) 
7 x 



10 g cm . This represents a small fraction (between 
0.7 and 4%) of the total cocoon density, which is in the 
range of (2 — 10) x 10 -32 g cm -3 , for inclination angles 
70° < 8 < 80°. Therefore, a substantial fraction of the 
mass density in the relict cocoon is IGM material ingested 
into the cocoon at the hot spot, and in instabilities along 
the contact discontinuity. At the same time the density 
in the cocoon is still less than the density of the IGM, 
which we estimate to be ~ 6 x 10~ 31 g cm -3 . 

This calculation highlights the fact that entrainment 
of the external IGM into the outer radio lobes is indeed 
significant. Consequently, entrainment must play an im- 
portant role in determining the evolution of the restarted 
jets in the r e lict sy nchrotron cocoons, as suggested by 
iKaiser et al.l (|2000T ) and as required by the theoretical 
interpretation given in this section. 



6 CONCLUSIONS 

The radio continuum images presented herein of the 
restarting radio galaxy PKS B1545— 321 show the inter- 
action between the new jets and relict cocoon plasma 
with the highest sensitivity and resolution to date. Most 
remarkable is the detection of a straight 3.2-kpc diameter 
jet along the axis of the inner SE lobe, which has a di- 
ameter of 20 kpc. This jet is observed not as an emission 
feature, but as a decrement in total intensity, presum- 
ably owing to a Doppler dimming of its emission and/or 
a lack of accelerated particles prior to the jet plasma be- 
ing shocked at the terminus. 

We do not detect any evidence for bow shocks asso- 
ciated with the restarted jets external to the inner lobes 
in PKS B1545-321. The Clarke et al. simulations pre- 
dict that such features ought to be visible with a bright- 
ness comparable to the hotspots. The lack of such ex- 
pected features, combined with the neat rims of emission 
at the ends of the inner lobes, lead us to propose that the 
boundary of the inner double is itself the bow shock of 
the restarting jet. Specifically, we suggest that the inner 
double lobes represent a mixture of relict cocoon ma- 
terial and new jet plasma that has been reaccelerated 
at the bow shock and terminal jet shocks respectively. 
It is along these lines that we also predict that there 
might exist giant radio sources which show no evidence 
for spectral aging along their axis, and that these sources 
may have recently restarted and had their entire cocoons 
reenergized via passage through a bow shock. Indeed, the 
classic discrepancy between dynamical and spectral ages 
in radio sources might be due to restarting. 

There are a number of observational features that 
support our above interpretation of the inner lobes in 
PKS B1545— 321: (1) As mentioned already, we have de- 
tected and resolved relatively bright emission peaks at 
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the ends of the inner double together with rims of emis- 
sion along the boundary. This morphology is uncharac- 
teristic of classic hot spots which usually protrude be- 
yond the associated lobe. These rims appear to be consis- 
tent with the emission from bow shocks propagating into 
the background non-thermal plasma. (2) The magnetic 
field is aligned with the rims of the inner lobes consis- 
tent with the alignment of magnetic field in a bow shock. 
(3) Local minima in the polarization of the inner lobes 
are coincident with or adjacent to the brightest features. 
Again, this is uncharacteristic of normal hot spots. We 
propose that these features are the result of the superim- 
posed emission from moderate strength Mach disks and 
bow shocks, with mutually perpendicular magnetic fields, 
which we expect from supersonic jets terminating in rel- 
ativistic cocoon material that is light compared to the 
normal IGM. 

In addition, we have performed a number of detailed 
dynamical calculations, which together with the new ra- 
dio continuum images, support the view that the evolu- 
tion of the new jets in the Mpc-scale lobes is unusual. 
We summarise these as follows and also note that each 
point is consistent with the above interpretation of the 
inner lobes as representing both bow shocks and new jet 
material: 

(i) The emission peaks at the ends of the inner lobes 
have a substantially lower surface brightness compared to 
that observed in the hotspots of powerful radio galaxies. 
This is indicative of the jets terminating in a low density 
environment. 

(ii) Treating the arm length asymmetry in the in- 
ner lobes as the result of time retardation, we have de- 
duced mildly relativistic hot spot velocities in the range 
0.3 < /3hs < 0.6 for the inner lobes, consistent with the 
propagation of the jets in a low density environment, but 
not as low as expected from pristine jet material. 

(iii) By considering the hot spot dynamics we have de- 
duced a density in the lobes which is less than the likely 
IGM density but which also exceeds the likely density 
produced by the original jets. The Mach number of the 
advancing hot spot is at least 5 with greater values pos- 
sible if the minimum energy hot spot pressure is signifi- 
cantly exceeded. 

(iv) The hot spot Mach number inferred by associating 
the spreading rate of the SE inner lobe with the Mach 
angle is approximately 23. 

(v) Consideration of the timescales related to the ex- 
pansion of the relict hot spots and the propagation of 
the new hot spots provide a lower limit on the advance 
speed (/3hs > 0.27) which is consistent with the range of 
velocities quoted above 

The standout discrepancy in the above is between 
the two estimates of the Mach number and we have ad- 
vanced three possible explanations for this. The first is 
that the jet wanders slightly in its direction producing 
a broader working surface thereby decreasing the appar- 
ent Mach angle. The second is that the expansion of the 
inner lobes is further impeded by a density gradient to- 
ward the sides of the cocoon. The third is that the hot 
spot pressure significantly exceeds the minimum energy 
value. The first idea may imply that the velocity of ad- 



vance of the jet is not constant as we have supposed. It 
is possible, for example, that the jet evolves according to 
the self-similar velocity oc £~ 2,/ 5 law for a con s tant b ack- 
ground density in the mode ls of lBicknell et al.l (|l997t ) and 
iKaiser fc Alexander! (|l997T ). Against this is the apparent 
lack of self-similar spatial structure in the morphology of 
the inner lobes. The self-similar models imply that the ra- 
tio of the rate of lateral expansion to forward expansion 
is constant and the base of the lobes (nearest to the core) 
should therefore have expanded more than is apparent in 
the images. Moreover, unpublished simulations by us of 
restarting jets show that the new jet is relatively unde- 
tected by the much lower density background medium 
of the pre-existing lobes. Nevertheless, it is possible that 
when adequately perturbed at the base such a jet may 
operate as a "dentist drill" but may not propagate as 
described by the self-similar models; this possibility de- 
serves attention, but detailed investigation is beyond the 
scope of this paper. 

There is support for the second idea that the jet 
may be affected by local density gradients, in the VLA 
and ATCA image in Figure 1. The extension of the trail- 
ing part of the NW inner lobe to the west indicates that 
the northern cocoon exerts a dynamical influence on the 
NW inner lobe. Also, the gradual curvature of the in- 
ner lobes points to a dynamical influence of the cocoon 
on the lobes. Further simulations of the interaction of 
a restarting jet with a pre-existing lobe should provide 
some valuable information on these issues. 

Notwithstanding this discrepancy between the esti- 
mates of the Mach numbers of the hot spots, there are 
several other features of the observations and our anal- 
ysis which are consistent with previous ideas on double- 
double radio sources. The deduction of mildly relativis- 
tic hot spot vel ocities is qualitatively consistent w ith 
the early models (IClarke fc Burnslll99ll ; IClarkelll997i ) al- 
though those models would predict velocities closer to the 
speed of light and fai nter inner lob e emis sion. Therefore 
we find the notion by IKaiser et al.l l|2000h , that entrap- 
ment of the surrounding IGM could decrease the sound 
speed in the cocoon, slow down the speed of advance 
of the inner doubles and make the inner lobe emission 
brighter, to be an appealing one. 

The interruption to jet activity in PKS B1545— 321 
has been brief. We have shown in <j5.6l that given our esti- 
mates of cocoon density, it is likely that the interruption 
to jet activity has been no more than a few percent of the 
age of the whole source. This figure could be less than a 
percent for a source lifetime ~ 2 x 10 s yr. How likely is 
it then that we would view this source at this stage in 
its lifetime? Effectively the time available for viewing a 
restarting jet is given by the travel time of the hot spot 
from the core to the extremity of the source « L/c/3hs in 
the notation of § 15.61 For the Northern lobe this is of order 
(3-4) x 10 6 yr and this could be (1.5-10)% of the source 
age, which we have estimated as being ~ (0.3 — 2) x 10 s yr. 
It is unsurprising therefore that double-double sources 
are relatively rare. Nevertheless, it would be of interest 
to examine the statistics of an unbiased sample of Mpc- 
scale radio galaxies to see if these sorts of fractional time 
estimates are appropriate. 

Finally, we note that we have shown that the host 
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Table 1. Journal of VLA observations 



Array 


Frequency 


Date 


Duration 




(MHz) 




Chr) 


BnA 


1384 


2002 May 31 


6 


BnA 


4910 


2002 Jun 2 


6 


CnB 


4910 


2002 Sep 27 


6 



galaxy is located at the boundary of a large scale filamen- 
tary structure, and shows blue patches in color distribu- 
tion indicative of a recent merger, presumably a minor 
merger, and it may be that the merger event triggered the 
Mpc-scale radio galaxy. Dating the blue stellar popula- 
tions that we have discovered in the host galaxy would be 
useful in estimating the time-delay between the merger 
and its influence on the emergence of the radio source. 
However, we should emphasize that the timescales we 
have estimated for the source switching off and then on 
again are not comparable to merger timescales ~ 10 8-9 yr 
and that one has to investigate other processes that would 
lead to such comparatively small time scales. 



APPENDIX A: IMAGING COMBINED VLA 
AND ATCA DATASETS 

Our VLA and ATCA 22 cm visibilities were combined 
and imaged using mosaic imaging techniques. The ob- 
servations made by the individual telescopes were not 
mosaic observations, in which visibility data correspond- 
ing to multiple antenna pointings are acquired. Rather, 
the VLA and ATCA observations were both single point- 
ing observations, and both observations used the same 
antenna pointing position. Nevertheless, it is necessary 
to treat the data specially while imaging because the 
antenna sizes differ in the two telescopes and, conse- 
quently the primary beams are different. Therefore, the 
sky source represented in the visibility data differ in the 
two telescopes, and sources offset from the pointing po- 
sition have different primary beam attenuations. 

The combined VLA and ATCA 22 cm image of 
PKS B1545— 321 shown in Figure [1] was computed in 
three separate steps, using the mosaic imaging, decon- 
volution and restoration tasks available in the MIRIAD 
package. Firstly, a linear mosaic of the 22 cm VLA and 
ATCA visibilities was performed using task INVERT. Af- 
ter converting the ATCA linear and VLA circular polar- 
ization measurements into Stokes parameters, the routine 
images the two pointings (visibility datasets) separately 
and then combines the resulting dirty images in a lin- 
ear mosaic process. Pixels in the individual images are 
weighted to correct for the primary beam attenuation, 
as well as to minimize the noise in the combined image. 
Specifically, the resultant dirty image that is output by 
INVERT, I(t,m), is given by 



I(£,m) = W(£,m) 



)Ii(£,m)/a? 



(Al) 



^Pfil-lum-mi)/*? 

(|Sault et alj|l996h where the summation, i, is over the 
pointing centers (fj,m,). The parameter Ii(£,m) refers 
to the image formed from the i'th pointing and Pi(£,m) 
is the primary beam pattern for the corresponding tele- 
scope. W(£,m) is a weighting factor, which ensures that 
the noise at each mosaic pixel does not exceed a particular 
threshold value. In the case considered here, the pointing 
centers refer to the individual ATCA and VLA observa- 
tions and the weighting function, W(£, m) = 1, within the 
region of the mosaic containing PKS B1545— 321. There- 
fore, the attenuation due to the VLA and ATCA primary 
beams is totally corrected for over this area. 

In addition to creating the combined dirty image, 
INVERT outputs a data cube containing the dirty sythe- 
sized beam patterns and primary beam models for each 
pointing. This enables the true point spread function, 
which is essentially a linear mosaic of the individual syn- 
thesized beam patterns, to be computed at any position 
in the image during the deconvolution stage. 

In the second step, the MIRIAD routine MOSMEM 
(|Sault et alj|l996h was used to perform a joint maximum 
entropy deconvolution of the 22 cm mosaiced dirty im- 
age. For each pointing, the routine multiplies the relevant 
primary beam with a prospective global model of the sky 
and then convolves the product with the dirty beam for 
that pointing. The dirty images are then linearly mo- 
saiced together, and subtracted from the combined dirty 
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RA (J2000) 

Figure 1. The combined VLA and ATCA 22 cm image of PKS 1545-321 made with a beam FWHM 3'.'6 X 2" 9 at a PA of 54°. 
Contours are at (-1, 1, 2, 3, 4, 6, 8, 12, and 16) X 200 fijy beam -1 . Grey scales are shown in the range 0.2-4.0 mjy beam - using 
a linear scale. The rms noise in the image is 40/xJy beam -1 . The half-power size of the synthesized beam is displayed in a box in 
the bottom right-hand corner. This image, as well as all others displayed herein, has been corrected for the attenuation due to the 
primary beam. 



image, to form the residuals. Through this process, the 
position-variant point spread function of the linear mo- 
saic is accounted for. 

In the final step, the global sky model and deconvo- 
lution residuals, as determined by MOSMEM, were input 
to the MIRIAD task RESTOR to compute a clean im- 
age. RESTOR convolves the global model with a constant 



Gaussian CLEAN beam (a beam that is not a function 
of position) and then folds the residuals into the result. 
The actual point spread function was calculated at var- 
ious positions in the linear mosaic; these were found to 
be fairly constant with parameters similar to that of the 
synthesized VLA beam. Therefore, we let RESTOR fol- 
low the default approach, which is to convolve the model 
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RA (J2000) 

RA (J2000) 

Figure 2. The southern and northern inner sources of PKS 1545—321 at 6 cm wavelength made with a beam FWHM l'/3 X l'/3. 
Contours are at (-1.0, 1.0, 2.0, 3.0, 4.0, 5.0, 5.5, 6.0, 6.8, 8.0, 9.0, 10.0, 11.0 and 12.0) X 65 £tjy beam' 1 . Grey scales arc shown 
in the range 0.065-1.0 mjy beam -1 with a linear scale. The rms noise in the image is 20 /ijy beam -1 . The half-power size of the 
synthesized beam is shown in the bottom right-hand corner. 




15 h 4B m 56 s 55 s 5 55 s 54 B 5 54 s 

RA (JS000) 

RA (J2000) 

Figure 3. The southern and northern inner sources of PKS 1545—321 at 22 cm wavelength made with a beam FWHM 2'.'6 X 2'.'6. 
Contours are at (-1.0,1.0,2.0,3.0,3.8,4.2,5.0,6.0,8.0,10.0,12.0) X 400 /xJy beam" 1 . Grey scales arc shown in the range 0.2-5.0 mjy 
beam -1 with a linear scale. The rms noise in the image is 80 fj,3y beam -1 . The half-power size of the synthesized beam is shown 
in the bottom right-hand corner. The right-hand panel also shows a slice profile across the NW inner lobe taken along the line 
indicated by the arrow. 
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Figure 4. Polarization over the southern and northern inner double. The images have been made with a beam FWHM 3"5 
X 3'.'5. The 6 cm polarized intensity is shown using grey scales. Bars show the B-field orientation with lengths proportional to 
the fractional polarization at 6 cm; the vectors shown in the bottom right-hand corners correspond to 100% polarization; the 
orientations of the B-vectors have been corrected for the line-of-sight Faraday rotation. Contours of the 6 cm total intensity are 
overlaid; contour levels are (-1, 1, 2, 3, 4, 6, 8, and 12) X 65 /^Jy beam -1 . 



with a CLEAN beam whose parameters are determined 
by fitting to the synthesized beam of the first pointing; 
this was the VLA pointing in our case. 



APPENDIX B: DISSIPATION TIMESCALE 
FOR A HOT SPOT 

If the beams from the central engine stop abruptly, and 
the last jet material passes through the hot spots at the 
ends of radio lobes, the hot spots will dissipate into the 
surrounding lobe. In this section we estimate the dissipa- 
tion timescale assuming that the hot spots are spherically 
symmetric and the expansion is adiabatic. 

Let the initial and final radii, pressures and rest mass 
particle densities within the hot spots be Ro and R, po, 
Phs, po and p, respectively. The thermal particle density 
in the external lobe is denoted by p ox t. Assuming that 
the expansion of the hot spot, with speed v cxp , is ram 
pressure confined as a result of the entrained thermal gas 
in the external lobe, 



Phs 
pcxt 



1/2 



(Bl) 



Since pR 3 is a constant in the expansion, which we as- 
sume is adiabatic, 



(B2) 



Phs _ I P 

Po {.Po, 

where 7 = 4/3 is the adiabatic index of the relativistic 
gas within the expanding hot spot. It follows that 

tH^r <b3 > 

The expansion velocity is then 

dR _ {PhA i/2 (po_Y /a (]L 



dt 



K Pcxt 



\Pext / \Ro 



(B4) 



Denoting R/Ro by y, 



1/2 



dy _ ( Po 
dt \p C xtRl 

This equation has the solution, 



V = 



1 + 3 



( Po 



1/2 



1/3 



(B5) 



(B6) 



which may be rearranged to give an expression for the 
expansion time t: 



t = 



-3/4 



pc^tRo 



po 



(B7) 



APPENDIX C: ESTIMATE OF JET ENERGY 
FLUX, FOR A RESTARTED JET, FROM HOT 
SPOT PARAMETERS 

In this section we use the momentum balance at the end 
of the hot spot to express the jet energy flux in terms of 
the observable parameters of the hot spot and the density 
of the cocoon. 

As in the body of this paper, we let /3j e t be the ra- 
tio of jet speed to the speed of light as measured in the 
rest frame of the host galaxy, Ljct is the corresponding 
Lorentz factor, Aj et is the cross- sectional area of the jet. 
Phs, Ths and Ahs are the corresponding parameters for 
the hot spot. p c is the density of the cocoon surrounding 
the restarted jet. We also allow for the fact that near the 
working surface the momentum flux may be spread over 
a wider area than the instantaneous jet cross-section ei- 
ther as a result of the jet direction changing or the flow 
dynamics in the vicinity of the working surface. Hence the 
cross-sectional area of the hot spot A^b ^ Aiet- 

We denote the jet speed and corresponding Lorentz 
factor in the frame of the contact discontinuity, or hot 
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Figure 5. Distribution of spectral index over the inner source computed from VLA images at 6 and 22 cm with beams of FWHM 
3". 5 X 3". 5. Grey scales are shown in the range of 0.1 to —1.0. Contours of the 6 cm total intensity are overlaid; contour levels are 
(-1, 1, 2, 3, 4, 6, 8, and 12) X 160 /xjy beam" 1 . 



spot, by additional subscripts: 'Its'. In the frame of the 
hot spot, momentum balance gives: 

rj 2 c t,h s /3fot,h s (Pjotc 2 + 4p jc t)Aj ct = p c c 2 pl B A hs . (CI) 

Lorentz transformation yields a relationship between pa- 
rameters in the frame of the host galaxy to that in the 
frame of the hot spot: 

rfet.hs^jet.hs = Thal^t (Act — (3 hs ) . (C2) 

and the momentum balance equation reads 

4rLrj 2 ot f>jct(l + Xjct)(/3jot - Ai S ) 2 Ai<=t = PcC /3 2 s A hs (C3) 

where Xjet = Pj ct c 2 /4pj ot . 

The energy flux in the relativistic beam is given by: 



F E = 4 C rf ct /3 jctPjct ( 1 + ^ - Xjat ) A jet . 



(C4) 



Dividing this equation by equation (|C3|) gives the follow- 
ing expression for the energy flux: 



F E = 



1 + - 



"Xi<=t 



1 + Xjc 



x/9 C c i3 jet r^ 



■ (/3 hB //3 jet ) 2 
[1 " (/Wftet)] 



:A hs .(C5) 



For even modest Lorentz factors for the jet, the leading 
factor in square brackets is close to unity, irrespective 
of the value of Xjct- Additionally, /3j et is close to unity. 
Hence this expression may be used to estimate the jet 
energy flux in terms of the hot spot speed, hot spot area 
and cocoon density. 
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Figure 7. B-band AAT WFI image of the host galaxy. 
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Figure 9. Residual B-band AAT WFI image of the host galaxy after subtracting the model galaxy. The central 2 arcseconds are 
not included in the model. 
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Figure 10. Panels (a) through (e) show the smoothed 6dF galaxy density (per sq. degree) in the vicinity of the host galaxy of 
PKS B1545— 321. Panels are 20 Mpc deep in redshift space and overlap by 10 Mpc. Contour levels are at: 0.1, 0.2, 0.4 and 0.8. The 
open circles mark the locations of 6dF galaxies in the corresponding redshift slices. The position of the host galaxy is indicated 
with a star symbol. Panel (f) shows the redshift completeness of the 6dFGS. 
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Figure 12. Hot spot advance velocity (3^ as a function of inclination angle 8, implied by the arm-length asymmetry ratio of 
D = 1.22 ± 0.02. The dashed curves show the lower and upper limits respectively implied by the uncertainty in the arm-length 
asymmetry measurement. 
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Figure 13. Density in the outer lobes p c as a function of inclination angle 
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Figure 14. The external Mach number of the hot spot M^ s as a function of inclination angle 8 for p^a = Pmin (solid line) and 
lOpmin (dot-dashed line). 
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Figure 15. Plotted widths of profiles of the SE trail. 
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Figure 16. The energy flux of the restarting jets Fj; as a function of inclination angle 6, computed from equation 1121 (solid line) 
and equation 1141 (dot-dashed line). 
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